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I. INTRODUCTION 



Studies on microscopic models for vehicle traffic provided a good point of view on the 
phase transition from the free traffic flow to the congested one. Related self-driven many- 
particle systems have attracted considerable interests not only from engineers but also from 
physicists jl, 2|. Among such models, the optimal velocity (OV) model which is a 
car-following model describing an adaptation to the optimal velocity that depends on the 
headway between two neighboring vehicles, is well-known for its successful rationalization 
of "phantom traffic jams" in the high-density regime. 

Whereas the OV model consists of ordinary differential equations (ODE), cellular au- 
tomata (CA) such as the Nagel-Schreckenberg model 41, the elementary CA of Rule 184 

n n n 

(ECA184) [5|, the Fukui-Ishibashi (FI) model |6| and the slow-to-start (s2s) model [7] are 
extensively used in analyses of traffic flow. A way toward amalgamation of the ODE-type 
and the CA-type models was opened by the discovery of the discrete OV (dOV) model [8] 
that provides an ultradiscretization 9] of the OV model. The resultant ultradiscrete OV 
(uOV) model includes both the ECA184 and the FI model as its special cases. We should 

n 

note here that another ultradiscretization |10| of the OV model is developed out of the 
ultradiscretization of the mKdV equation ll|, which also has an application to the traffic 

flow bj. 

n 

Inspiration brought by the ultradiscretization of the OV model in ref. [8J leads us to 
an ultradiscretizable hybrid of the OV and the s2s models, which was named the optimal 

n 

velocity model with the slow-to-start effect, or shortly, the s2s-OV model |13l |. 
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where to and Xk(t) are the monitoring period and the position of the k-th car at the time t, 
respectively. The interval between the cars k and k+1 is denoted by Axk(t) = Xk+i(t)—Xk{i). 
Sensitivity to the interval is controlled by 5x > and xq means the length of the road 
occupied by a vehicle. In the limit in — > 0, the s2s-OV model reduces to the Newell 

n 

model [14J , which is a car- following model describing an retarded adaptation to the optimal 
velocity determined by the headway in the past. The ultradiscrete limit of the s2s-OV model 
includes a traffic-flow model of the CA-type, 
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where the position of the k-th car x^, k = 1, 2, • • • , K, and the headway between the cars k 
and k + 1 at the n-th discrete time Ax^ := x% +1 — x% — xq > are integers. We should note 
that the headway Ax£ and the interval Ax\ := x"^ +1 — x\ = Ax\ + Xq between the cars k 
and k + 1 are used in different meanings in this paper. If all the initial headways Ax° k are 
non- negative, then car crash is prohibited at any time. Positive integers ^o, xq and 5t means 
the speed limit, the size of a cell and the discrete time step, respectively. When we consider 
the time-evolution of the vehicles in the following, the latter two parameters, Xq and 5t, are 
fixed at one. Then eq. ([2]) reduces to 

"0 



k — 



£ + min(min(*— - x\~ " -l),v j. (3) 



We call a positive integer uq the monitoring period since the minimum of the headway 
Ax^~ n := SfeT" — x k~ n — %o of n + I discrete times, n — n' = n — n , n — n + 1, • • • , n, 
is involved in the time evolution ([2]). This CA describes many cars running on a single lane 
highway in one direction, which is driven by cautious drivers requiring enough headway to 



erate their cars. Since the above CA (T5]) is 
incorporating the s2s effect in the sense 



go on at least for uq time steps before they acce 
an extension of the CA case of the uOV model 
of delay in drivers' response, we call the CA ([2]) the CA variant of the OV model with the 
s2s effect, or the s2s-OVCA for short. 

Figure [1] gives typical examples of the spatio-temporal pattern showing jams and the fun- 



damental diagram of the s2s-OVCA 13|. The spatio-temporal pattern shows the trajectories 
of vehicles. Though irregular patterns are observed at the initial stage of time evolution, 
trajectories become regular in the long run, which consist of several jam clusters moving 
backward at the same velocity and vehicles running at the top speed. The fundamental 
diagram, or the flow-density relation, in other words, gives the relation between the vehicle 
flow 

1 JL T n+1 _ r n 
q ._ t \ ^ \^ ±k ±k_ 

(n n . + 1 ) L Z^Z^ st ' 

v ■> ' k=l n=ni 

which is equivalent to the total momentum of vehicles per unit length averaged over the 
period n« < n < rif when the traffic is expected to be in equilibrium, and the vehicle density 
p := where L is the total number of the cells, or the length of the circuit, on which the 
periodic boundary condition is imposed. The fundamental diagram has the inverse-A shape 
with several phases of traffic, namely, free, jam as well as metastable slow traffics, which 



captures the characteristic of observed flow-density relations 



2J. A comment here might 
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FIG. 1. The spatio-temporal pattern (left) and the fundamental diagram (right) of the s2s- 
OVCA [2]. The number of the cells L is fixed at L = 100 and the periodic boundary condition is 
imposed. The number of the cars in the above spatio-temporal pattern is K = 30. The maximum 
velocity vq and the monitoring period no are vq = 3 and no = 2. The flows Q are computed by 
averaging over the time period 800 < n < 1000, in which the traffic is expected to be in equilibrium. 



be in order. A CA-type model that is different from the s2s-OVCA showing a fundamental 
diagram with several branches was reported. 15] But all of its branches except for the free 



line have negative inclinations. This feature is also different from that of the fundamental 
diagram of the s2s-OVCA. As we can see in Fig. (TJ the fundamental diagram consists of 
straight branches, whose number, four in this case, is the same as that of all possible integral 
velocities, v — 0,1,2 and 3(= Vq). The main result of this paper is to present a set of exact 
solutions of the s2s-OVCA, which explains the piecewise linear fundamental diagram of the 
s2s-OVCA. 



The outline is as follows. In ^TTJ the fundamental diagram given in Fig. [T] will be examined. 
An empirical formula for the fundamental diagram of the s2s-OVCA will be introduced. 
A set of exact solutions of the s2s-OVCA that reproduces the empirical formula will be 
presented in §1111 and §IV1 Concluding remarks are given in the final section. Derivation of 
the s2s-OVCA (J2]) from the s2s-OV model ([TJ is briefly summarized in an appendix. 
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II. NUMERICAL OBSERVATION OF THE FUNDAMENTAL DIAGRAM 



By numerical experiments 13| , the fundamental diagram of the s2s-0VCA was observed 
to be piecewise linear and possesses v branches irrespective of the parameter n . Let us 
have a close look at the fundamental diagram in Figs. [T] and [TT] and examine these features 
more in detail. 
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FIG. 2. (Color online) An empirical formula for the fundamental diagram of the s2s-OVCA. All the 
dots on the diagram were numerically computed [13j for the case in which the maximum velocity 
vq and the monitoring period no are v$ = 3 and uq = 2. Each solid line corresponds to the formula 
in the same color or in the same grayscale. Each branch is labeled with a number that shows the 
integral velocity v corresponding to the branch. The end points of the branches are determined 
by the maximum density for each integral velocity. The branch points of the slow or the jamming 
lines, which are encircled with thin small circles, give the lowest densities for these lines. 

As shown in Fig. [Til each branch corresponds to an integral velocity that is less than or 
equal to the maximum velocity vq. All the end points of the branches are observed to be on 
the line 

Q5t + px = 1. (4) 

This is because the density at the end point is determined by the maximum density p max (^) 
that allows the velocity of the slowest car to be v. Because of the time evolution of the 
s2s-OVCA d2J), a car running at a velocity v must have an interval that is larger than or 
equal to v5t + Xq. Imagining a traffic such that all the cars flow at the same velocity v with 
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the minimum interval vSt+xo, we conclude that the maximum density Pmax(w) is determined 
by 

Pma X (f ) = -f— , (5) 

vdt + x 

since each car occupies vdt + x$ cells. The flow corresponding to the traffic is given by 

Q(Pwax(v)) = Vp max (v). (6) 

Equations (jSJ) and (jSJ) simultaneously holds on the end points and they yield Q (p max (v)) 5t + 
p max (f)x = 1. Thus all the end points must satisfy the relation (J3J). 

The free line is a branch that agrees with the straight line whose inclination equals to 
the maximum velocity v Q , 

Q = v p, (7) 

since all the cars on the free line flow at the maximum speed vq. Any other branches branch 
off from the free line. By observation, the density of the branching point of the branch 
corresponding to a velocity v reads 

Pmm(v) = —, \~rT j TT, • ( 8 ) 

n>o{v - v jot + v ot + x 
This observation is explained as follows. Suppose one car, say the car k, runs at the velocity 
v and all the other K — 1 cars run at the maximum velocity Vq. At the moment the k-th 
car slows down to v, the interval between the cars k and k + 1 is vSt + Xq. Since it takes 
at least hq + 1 time steps for the car k to speed up to vq, the interval between the cars k 
and k + 1 expands up to H = (no + l)(fo — v)5t + vSt + xq = l/pmm{v) > v o$t by the time 
the k-th car speeds up to Vq. If all the cars can obtain the interval H, slow cars running 
at the velocity v should vanish in the end. Thus the density at the branching point of the 
branch corresponding to the velocity vSt is given by p m i n (v) = 1/H, which is nothing but 
the minimum density of the branch. 

Suppose the slow branches form straight lines. Then these lines must be the segments 
connecting the branching point (pmi n (v), UoPmin(f )) and the end point (p max _(v),vp maK (v)) . 
The expression for these line segments is obtained by eqs. (jSJ), © and (jSJ) as 

_ n vSt - x 1 , v 

V " (no + l)5t P+ {no + lW [ ' 

which agrees fairly well with the branches for the slow traffics consisting of the numerically 
obtained dots in Fig. [TXJ Coincidence of the Q- interception of slow branches ^ no ^ St , which 
also appears in Fig. [TXJ, is also shown in the above empirical formula. 



The aim of the following two sections is to explain the empirical formula for the funda- 
mental diagram fl9]) based on exact solutions of the s2s-OVCA. 



III. SOLUTIONS WITH A SINGLE SLOW CLUSTER 

After transient flow goes on evolving irregularly for a while, the traffic reaches "equilib- 
rium" in the long run where slow-traffic clusters of seemingly invariant sizes move regularly, 
as one can see in the spatio-temporal pattern in Fig. [TJ Here we allow equilibrium to change 
sizes of slow or fast clusters periodically but prohibit it from increasing or decreasing their 
sizes monotonically. We shall investigate such equilibrium flows and their flow-density re- 
lations. Limiting ourselves to simple cases, we begin with equilibrium flows with, at most, 
a single slow cluster. 

To begin with, we shall see a specific spatio-temporal pattern starting from the following 
initial configuration, 

: 1u2u3u4uuuuuuu5uuulj- (10) 

Note that the number at the leftmost shows the time. The digits and the blank symbols 
u in the above configuration mean the indices of the cars and the empty cells, respectively. 
We set the monitoring period no at 2. The velocity of the cars 4 and 5 is 3, which is the 
top velocity vq of this case. The velocity of the cars 1, 2 and 3, whose headways are 1, is 1. 
Then in this case, the headways of all the cars in the past have nothing to do with the time 
evolution starting from the above configuration ([TO]) . The periodic boundary condition is 
assumed so the length of the circuit is 19 in this case. 

Using eq. (j3J), which is equivalent to the s2s-OVCA (j2J) with xq — St — 1, we can confirm 
that the following spatio-temporal pattern evolves out of the above initial configuration (JTOl) . 

: 1u2u3u4uljuuuuu5uuuu 

(11) 

2 : 5|jlu2u3uuuuu4uuuuuu 

3 • u5u1u2u3uuuuuuu4uulj- 

Moving all the cells of the initial configuration one cell rightward as well as neglecting the 
difference of the car indices, we get the configuration at the time 3. In this sense, the initial 
configuration ffTOl) gives a periodic motion of the vehicles and hence an equilibrium flow of 
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the s2s-0VCA (J2J). Roughly speaking, we shall construct equilibrium flows of the above 
type in the following. 

The most simple equilibrium flow is the free traffic in which all the cars always run at 
the top velocity Vq. In other words, the free traffic is a traffic flow with no slow cluster. 
The headway of a car running at the top velocity must be larger than or equal to v^St and 
so is each headway in the free traffic. Since all the cars share the same velocity in the free 
traffic, all the headways are invariant. Thus the free traffic is in equilibrium. One concludes 
that the flow-density relation for the free traffic is given by the free line (J7|) because of its 
definition given above. 

A snapshot of an equilibrium flow with only a single slow cluster is given by a configuration 
of vehicles shown below, 

,\ [01 ^[0l/ x [0], fe+1 x [n +l] f-hno+1] K [m+l]-lx 1 

(«)-*■ ■ ■ ■ (v)-t{v)~*(v ) -> ■■■{v ) -»■ {v ) ~» I (v), (12) 



L 

where m = 0, 1, ■ • • ,hq and lxo/dt = 0, 1 • • ■ ,Vq — v — 1. The car n, n — 1, 2, • • • , K, which 

n 

will move at a velocity v, v = 0, 1, • • • ,vq, in the next turn, is denoted by (v). Thus the cars 
1 to A; form a single slow cluster in the above configuration and the cars k + 1 to K form 
a single fast cluster. Thus the integer k = 0, 1, • • • , K gives the number of cars in the slow 
cluster. We shall not deal with the case k = since its corresponding configuration gives a 
free traffic without slow clusters. Thus k is limited to k — 1, 2, • • • , K in the following. 

An arrow A,/i = 0,l,2,---,is the headway whose number is h, or h empty cells between 
the neighboring cars. The symbol [n] that appears on the arrows means [n] = v5t + ti(vq — 
v)St. Two headways in front of the cars k and K, which are the tops of the slow and the fast 
clusters respectively, are denoted by wavy arrows ~^ in order to stress that they are going to 
change by (vq — v)5t at each step. Dots between a pair of cars with the same velocity v and 
headway [h], (v)—f ■ ■ ■ (v)—t, mean that all the cars and the headways in between also have 
the same velocity and headway, (v)—t, n\ < v n < ni. The vertical bars mark the initial cell 
of the car 1. The car 1 and the vertical bar, which appear again on the rightmost, show 
the periodicity of the circuit and they will be omitted in the following. Note that the initial 
configuration (flTJ|) is given by the snapshot (|T2|) with t> = 3, v = 1, n = 2, K = 5, k = 3 
and I = 3, 

1 2 3 4 5 

: l u 2 u 3 u 4 uuuuuuu 5uuuu = | (1)^(1)^(1)^(3)4(3)- 
where 1 = [0], 7 = [3], 4 = [3] - 3 and x = 5t = I. 
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The circuit length L, which is shown under the brace in the above configuration, is then 
given by 

L = L Q — (no — tti)(vq — v)St — lx Q , (13) 

where L := (K — k)(n + l)(v — v)St + K(x$ + v St). With use of Pmin(^) given in eq. (jBJ), 
the formula for the circuit length flT3|) is cast into 

K 

k(n + l)(t>o — v)5t = — — L — (n — m)(t> — v)5t — Ixo, (14) 

Pining J 

where l/p m m{v) is an integer. Since the parameter k gives the number of cars in a slow 
cluster, it must be non-negative. Equation (fT4"l) gives the condition for k to be non-negative 

as 

K 

< — — L — (n — m)(vo — v)5t — lx$. 

Pmm{V) 

Thus we obtain the necessary condition for k to be non-negative as 

K 

P ■= J- > Pmin(^), (15) 

since ( n o- m )( v o- v ) St + lx o > q rj ng j ower i3 0un( j Q f p j n ^ ne above inequality coincides with the 
numerically obtained minimum density (JSJ for the branch corresponding to a velocity v. 

We should confirm here that there always exist non-negative integers k, I and m satisfying 
eq. f lT3|) for arbitrarily given L, K,Uq, Vq and v that meet the condition (1151) . The r.h.s. of 
eq. ( IT41) must be divided by (no + l)(t>o — v)5t. Since the third and fourth terms cover all 
the residue classes modulo (no + l)(v — v)St, there always exist some non-negative integers 
k, I and m satisfying the above equation for arbitrary L, K, uq, vq and v. 

The time evolution out of the initial configuration ( !T2|) gives an equilibrium flow with 
a single slow cluster. Let us confirm it by checking each step of the time evolution. We 
exclude the case k = K and shall deal with it separately. Watching the initial state again, 

: | (v)-J. • • • {v)-l{v)~>{vo) -> ■■•{v ) ->■ (vo) , (16) 

where at the leftmost shows the time, we notice that the cars 1 to k run at a slow velocity 
v and the cars k + 1 to K run at the top velocity Vq. Thus only two headways denoted 
by wavy arrows for clarity, namely those of the cars k and K that are the tops of the slow 
cluster and the fast cluster respectively, are going to change by (vq — v)St at each step. The 
configuration of the next time step is then given by 

1 [l]f +1 \K+l] f-\[n +l\ K [m\-lx 

1 : | =>{v)-$- ■ ■ ■ (tV-K^-^fo) -»■ ■■■{v ) -> (t> ) ~> • 



The subscript 1 of the car k shows that it is the first turn for the car k to have a headway 
which is long enough, [1] = vo5t, to run at the top speed v$. Similarly, a subscript n of the 

k 

k-th car, (v) n , v < vq, whose headway is longer than or equal to v 5t, means that it is the 
n-th turn for the car to have its headway long enough to run at the maximum velocity Vq. 
This symbol is going to appear shortly and to be used in the same meaning in the following. 
The thick arrow at the leftmost shows the displacement of the car 1 from its initial cell. 
Until the car K decelerates at the m-th time step, the cars 1 to k run at a slow velocity v 
and all the others run at the top speed vq. Thus the configuration at the n-th time step, 
1 < n < m — 1, is expressed by 

.nvSt.KlO] /"vPlA N, fc+ \N+1] f ~\ [n +l] , K . [m+l-n]-lx 

n : I =4> {v)-> ■ ■ ■ {v)->{v) n *^{vo) ->■ • • • {v ) ->■ (u ) ~» ■ (17) 

At the m-th time step, the headway of the car K shortens to [1] — Ixq = vo5t — Ixq and its 
velocity immediately slows down to v — lx /5t. Thus there additionally appears another 
varying headway denoted by a wavy arrow in front of the car K — 1. The configuration at 
the m-th time step is then given by 

m:| =4> (u)-^ • • • (u)-^ {v)^{v ) -> • • • (u ) ~* {v Q - lx a / 6t) ~* . (18) 

At this moment, velocities of the cars K — l, K and 1 are respectively v , v — lx /5t and v, as 
well as the headways of the cars K — l and K irregularly shortens by Ixq and (i>o — v)5t — Ixq- 
We classify the car K with the intermediate velocity vq — Ixo/St into the fast cluster. The 
slow cluster with the velocity v consists of cars with only a single slow velocity v. Thus the 
car K has been the top of the fast cluster by this turn. But in the next turn, the car K also 
slows down to v and becomes the bottom of the slow cluster. Thus the configuration at the 
(m + l)-th time step is 

(m+T)v8t 1 [01 ^jOl fc [m+1] fe+1 [n +l] f ~?[no+l]f -J[no+l]-teo A M 

Until the car k accelerates at the (n + l)-th time step, k + 1 cars, namely the cars 1 to A; as 
well as the car K, run at a slow velocity i> and all the others run at the top speed v . Thus 
the configuration at the n-th time step, m + 1 < n < n , is expressed by 

n:| =4> (v)^- • • • {v)-> (v) n ~^{v ) ->■ • • • (i»o) -> (^o) ~» ■ ( 19 ) 

At the (no + l)-th time step, the car k accelerates to the top velocity v since its headway 
has kept no less than v 5t for n time steps by the time. Thus the car k switches to the fast 
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cluster and the cars 1 to k — 1 as well as the car K now form the slow cluster. The place of 
the varying headway denoted by a wavy arrow also moves from the car k to the car k — 1. 
Thus the configuration at the (n + l)-th time step looks as 

(n +l)v8t 1 [0] k .-h0] . * [no+lh^lno+l] f "? [no+l]f ~y [m+l]-lx .\ [01 , . 

Since the vehicles are running on the circuit of the length L, we may put the car K behind 
the car 1. Then the configuration at the (no + l)-th time step 

. n v6t-x , K . [01 /. [01 [0] fc [no+1] f ~?[no+l]f -J[m+l]-te 

no + 1: | =>• (u>-^(t;>-J> • • • (u > -> • • • (uo> -> («b) ~» , (21) 

becomes the same as the initial configuration up to the rightward displacement of novSt — xo 
cells. Thus the motion of vehicles evolving out of the initial configuration (|16p is periodic 
with the period no + 1, which shows the corresponding traffic is in equilibrium. 
The configuration for the case k = K 



,nvSt,\[0\ , X , [01 

n : | =>- (i>)-4 • • • (v )—> 



gives a constant traffic flow in which all the cars run at the slow velocity v with the same 
headway [0] = v5t. Thus we confirm that the case k = K also gives an equilibrium traffic. 

We should note that the fast cluster in the case k = K — 1 temporally vanishes from the 
(m + l)-th time step to the n -th time step, 



: | (v)—i ■ ■ ■ (v)-i (v) £*{vq) 

iwrft/Vo] ^'xH^'x [n] , K x [m+l-n]-lx 



m+1]— IxQ, 



, rn^t/y [01 f \ 2 [0]f y 1 M/ f /r+vM-ko 



(22) 



n«5t , \ [01 5"\ 2 [0] ^ [ m + 1 l "'*o A [0] 



n : | '^(tj • • ■ (v) ^ (v) n l +> u (t>)-^-, m + 1 < n < n . 
But it appears appears again at the (n + l)-th time step 

(n + l)vSt 1 [01 ^[Olf-^Olf-hm+ll-teo,^, [0] 

n + 1 : | =>- (f)+-- - (v) — J- \v ) ~^(fo) ~^ {v}-^, 



which goes back to the shifted initial configuration 

| no^-zo/Vol/yM 5- 3 [0l^- 2 [0]f-l[m+l]-Zx o 



Thus such vanishment of the fast cluster as we have observed above does not break equilib- 
rium of the traffic. 
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Now let us compute the flow of this periodic motion averaged over the period no + 1. As 
we have already seen in the above time evolutions for the cases k — 1,2, • ■ ■ , K — 1, k and 
K — k cars respectively run at the velocities v and Vq in < n < m — 1. At the m-th time 
step, k and K — k — 1 cars run at the velocities v and Vq and remaining one car move at 
Vq — lxo/5t. In the following period, m + 1 < n < n , k + 1 and K — k — 1 cars move at the 
velocities v and vo, respectively. Thus the average flow Q is 

Q = j- ^ ^ L (m(kv + (K- k)v ) + (kv + (K - k - l)v + v - lx /8t) 
+ (n - m) ((k + l)v + (K-k - l)u )) 

TIqV 5t — Xq 1 

= (n + l)St 9 + (no + l)6t' 
which is the same as the empirical formula for the flow-density relation 0. For the case 
k = K, the flow Q and the density p are respectively given by 

Q = pv, p 



vSt + Xq 

A straightforward calculation shown below 

_ (no + l)v5t _ (n v5t - x ) + (vSt + x ) 

Q ~ P (n + l)5t ~ (n + l)5t 9 

nov5t — xo 1 
-p + 



(n + l)5t r (n + l)St 
proves that the above two quantities are related by the flow-density relation 0. 

The above formula is interpreted in a different and more intuitive manner. The config- 
urations ( TL6l) and (12 lj) shows the rightward displacement of the entire configuretion (fT6l) 
by n v5t — Xq cells in no + 1 time steps. The flow provided by this motion of the entire 
configuration gives the first term p of the above formula. Here we should remind 

ourselves of the leftward displacement of the car K by L cells, namely whole the circuit 
length, which is fictitiously introduced to make the shifted initial configuration ( [21]) from 
the real configuration ( 1201) at the (no + l)-th discrete time. In order to compensate the 
underestimation of the flow caused by this fictitious displacement, we have to add the flow 
corresponding to the rightward displacement of the car K by L cells in n + 1 time steps 

1 L 1 



(n + l)L 6t (n + l)<5t' 
which agrees with the second term of the above formula for the average flow. 
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IV. SOLUTIONS WITH MULTIPLE SLOW CLUSTERS 

As we have observed in the spatio-temporal pattern in Fig. [TJ several slow or jam clusters 
that share the same slow velocity coexist in equilibrium. Here we make solutions as such. 

Roughly speaking, we can make an equilibrium solution with multiple slow clusters by 
putting several snapshots of the equilibrium solutions in the previous section together. Let 
us see how this idea works via observation of a specific example. Thanks to the periodic 
boundary condition, we can move the leftmost empty cell of the configuration (iTTj) to the 
rightmost, 

ulu2u3uuu4uuuuuuu5u = | ^ 1u2u3ulju4uijuijuiju5lju 

Putting the above configuration and the configuration (flOl) 

Iu2u3u4uuuuuuu5uuuu = | (1)4(1)4(1)4(3)4(3)4 
together, we get a configuration with two slow clusters as 

1u2u3uuu4uuljuijulj5uli6lj7lj8lj9uuuuuuu0uuulj 

1 2 3 4 5 6 7 8 9 

= (1)4(1)4(1)^(3)4(2)^(1)4(1)4(1)4(3)4(3)4 (23) 

Using eq. d2J), we can confirm that the time evolution out of the above configuration gives a 
spatio-temporal patterns below, 

: Iu2u3uuu4uuuuuuu5uu6u7u8u9uuuuuuu0uuuu 

1 : u1lj2u3uuuuu4uuuljulj5lj6u7lj8uuu9uuuuuuu0u 

2 '. 0|j1|j2|j3|juuuuuu4ljuuu5u6u7u8ljuuuu9uuuuuu 

3 • u0ulu2uuu3uuuuuuu4uu5u6u7u8uuuuuuu9uuu- 

Moving all the cells of the initial configuration one cell rightward as well as neglecting the dif- 
ference of the car indices, we get the configuration at the time 3. Thus the configuration f[2"3"j) 
gives a periodic motion of the vehicles and hence an equilibrium flow of the s2s-OVCA 
In the following, we shall confirm that the above construction of equilibrium flows with 
multiple slow clusters works in general, too. 
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As building blocks of such solutions, we introduce symbols denoting slow and fast clusters 
in the snapshots of the equilibrium solutions f lT6|) . ffTT|) . f|T8|) and f fl9|) in the previous section, 

d i \ h hf nvdt < n < m 

n:\A(S}^{F}^, d=l ^ - - , (24) 

I (n — i)v5t — Xq m + 1 < n < n 
where (S) -w and (F) -w respectively denote the slow and the fast clusters with varying 

h s / 

headways ~^ whose lengths are h s j on their tops. For instance, the slow and fast clusters 
in the snapshot of the m-th time step ( Fl8|) are respectively given by 



h s A[o] k 



in 



(S) !$:=(v)^.-. (v)n (v) 

h f [n +l] K ~ 1 i,o l , 

(F) F J>:=(v ) h ■ ■ ■ (v o y'^'(v - lxo/St) 



"o+l], f , c ^[l]-lx 



where h s = [m] and hf = [1] — Ixq. Note that the car K with the intermediate velocity 
vq — Ixq 1 5t is classified into the fast cluster. We also note that we can always rearrange the 
order of the vehicles into the form of eq. fj24")) . thanks to the periodic boundary condition. 

An equilibrium solution with N slow clusters are made with the slow and fast clusters ( 1241) 
with the same monitoring period no, the slow velocity v and the top velocity vq. First, we 
prepare N pairs of slow and fast clusters in the configurations of N equilibrium solutions with 

h s - hf. 

only a single slow cluster, (Si) (F{) ~4 % = 1, 2, • • • ,N. We shall consider a configuration 
made by putting these pairs of slow and fast clusters in line, 

(S l ) ^ (F 1 ) "4 (S 2 ) ^ (F 2 ) ■ ■ ■ (S N ) ^ (F N ) h £ . (25) 

where the periodic boundary condition is imposed. Thus the total length of the circuit is 
the sum of the lengths of all the clusters. 

We observed that vanishment of a fast cluster ( BSD d° es n °t break equilibrium because 
the vanished fast cluster was revived by acceleration of the top car of the slow cluster behind. 
However, the top car of the fast cluster does not slow down if the slow cluster in front of 
it vanishes. Thus a vanishing slow cluster does break equilibrium. If slow down to the 
slow velocity v of the top car of any fast cluster (Fi) is no later than speed up to the top 
velocity Vq of the bottom car of the slow cluster (jSVfi) in front of the fast cluster (Fi), such 
vanishment of slow clusters is prevented. Thus we require the above "no vanishing slow 
cluster" (NVSC) rule to be fulfilled by any pairs of the fast and slow clusters, (Fi) and 
(S i+ i), in the configurations f l25|) . 
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We shall see below that the above configuration ( |25l) satisfying the NVSC rule gives 
an equilibrium solution with iV slow clusters. Speaking more specifically, we shall confirm 

h s . he, 

that the motion of vehicles belonging to a pair of slow and fast clusters (Si) ~4 (Fi) ~4, 
i — 1, 2, - • • , N, in the above configuration f!2^|) with multiple slow clusters is the same as 
that without any other pairs of slow and fast clusters we have investigated in detail in the 
previous section. 

Consider the time evolution of a fast cluster and its varying headway (Fi) i = 
1, 2, • • • , N. It is determined by the motion of the vehicle at the bottom of the slow cluster 
(S-i+i) in front of it, which is different from its original partner in the solution with 
a single slow-cluster (Si) behind it. However, thanks to the NVSC rule, the bottom of 

hs i+1 hf i 

(Si + i) keeps the slow velocity v until the top of (Fi) -4 slows down to v. Thus the 
motion of the bottom of (5V-t-i) is the same as that of the bottom of (Si) and so 

h f . 

is the motion of the fast cluster (Fi) ~4 until its top slows down to v. Let us now think 

h f . 

about the time evolution after the top of (Fi) --4 slows down to v. In the same manner we 
observed in the solution with only a single slow cluster in the previous section, the top of 

h f . 

(Fi) ^4 keeps the slow velocity v at least for no + 1 time steps. The time evolution of the 

hf. 

varying headway, which locates in front of the second car of (Fi) --4, and the motion of the 

hf. hf. 

remaining cars of (Fi) ~4 are now determined by the above motion of the top of (Fi) ~4. 

h f . 

Thus the time evolution of the fast cluster and its varying headway (Fi) — > remains the 
same as that with only a single slow cluster we observed in the previous section at least for 

hf. 

n + 1 time steps after the top of (Fi) ~4 slows down to v and at least for n + 1 time steps 
from the beginning. 

Let us turn ourselves now to the time evolution of a slow cluster and its varying headway 
(Si) i = 1,2, ••• ,N. The time evolution is determined by the motion of the vehicle 
at the bottom of the fast cluster in front of it, which is nothing but its original partner in 

hf. 

the solution with a single slow-cluster, (Fi) — >. We now know that the time evolution of 

h f . 

the fast cluster (Fi) ~* is the same as those with only a single slow cluster for uq + 1 time 

fts- 
steps. Thus the time evolution of the the slow cluster (Si) ~4 and its varying headway in 

the configuration with multiple slow clusters f )25l) is the same as that in the configuration 

with only a single slow cluster for no + 1 time steps. 

As we saw in the configuration with only a single slow cluster in the previous section, the 

time evolution of the slow and the fast clusters as well as their varying headways are periodic 
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with the period n Q + 1. Thus, also in the configuration with multiple slow clusters ( |25|) . the 
time evolution of all of the slow and the fast clusters as well as their varying headways 
are periodic with the period n + 1. To summarize, we have confirmed that the motion of 

h s - hj. 

vehicles belonging to a pair of slow and fast clusters (Si) (Fj) ~4, i — 1, 2, • - • , N, in the 
configuration (I25p with multiple slow clusters is the same as that without any other pairs 
of slow and fast clusters we have investigated in detail in the previous section and that it 
gives an equilibrium solution with N slow clusters. 

Let us examine the flow-density relation of the traffic evolving out of the configuration 
with iV slow clusters fl2"oT) . Let L iy K iy Qi and pi := Ki/Li be the circuit length, the number 
of vehicles, the mean flow averaged over one period, i.e. Uq + 1 time steps, and the density 

h a - hf. 

of vehicles for the configuration with only one slow cluster, (Si) ~4 (Fi) ~4, i = 1, 2, ■ • • ,N. 



In §IIIt we confirmed that Qi and pi satisfied the flow-density relation (J9J), 

_ n vSt - x 1 , 0(i \ 

Qi ~ (n + l)St Pi+ (n + l)5f W 

The circuit length L and the number of vehicles K of the configuration with iV slow clusters 
are then expressed by L := ^f =1 Li and K := Y2i=i-^i- The mean flow Q that is again 
averaged over the period as well as the vehicle density p of the configuration with iV slow 
clusters are expressed by 

^ N N 

Q = j LiQi = XiQi, 

i=l 

N 



L 

i=l 



K _ 



with Xi := Li/L. Using an identity ^2 i=1 Aj = 1 as well as the flow density relation for the 
configuration with a single slow cluster (T2"6"|) . one straightforwardly confirms 

n v5t -x \ fr, n ov8t - x 



(n + l)5t r (n + l)5t' 

1 1 



i=l 

A' 



(n + l)St {n + l)5t , 



i=l 

which concludes that Q and p also satisfy the flow-density relation (j9j) . Note that the relation 
can be confirmed in a more intuitive manner, as we have shown for the configuration with 
only a single slow cluster. 

16 



V. CONCLUDING REMARKS 



We have introduced a set of exact solutions of the s2s-0VCA and have explained the 
piecewise linear fundamental diagram of the model. The empirical formula for the flow- 
density relation of the s2s-0VCA ([9]) has been read out of the numerically obtained fun- 
damental diagram in Figs. [I] and [TTJ A set of exact periodic solutions that possesses only 
a single slow cluster has been presented, which has reproduced the formula for the flow- 
density relation of the s2s-OVCA (JUJ). With the use of these solutions, another set of exact 
periodic solutions with multiple slow clusters, which have been observed in the numerically 
obtained spatio-temporal pattern in Fig. [TJhas been introduced. And these solutions again 
have explained the flow-density relation of the s2s-OVCA (Q. 

We should comment about a couple of problems in our mind related to the present results. 
The solutions given in §1111 and §IVI provide a set of equilibrium flows which correspond all 
the points on the fundamental diagram of the s2s-OVCA. But we do not know if the set 
we have presented in §1111 and §IVI gives all possible equilibrium flows or not. There might 
be a solution giving an equilibrium flow that does not belong to the set. By numerical 
experiments [131 ]. we observed that time evolutions out of all the initial conditions we ver- 
ified eventually go to some equilibrium flows. But this observation still remains numerical 
observation and lacks proof. For the case of no = 1, these problems were investigated [3] 
and should be extended to all the cases of no's. As was shown in our previous paper jl3j . the 
s2s-OVCA 02]) transforms to the s2s-OV model ([T|) through the inverse ultradiscretization 
and the continuous limit. Note that the monitoring period no in the s2s-OVCA must be a 
free parameter so as to make the monitoring period t in the s2s-OV model a free parameter. 
It should be clarified if the characteristics of the equilibrium solution of the s2s-OVCA given 
in this paper remain in the solution of the s2s-OV model. 
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Appendix A 



We shall derive the s2s-0VCA (J2]) from the s2s-OV model (CQ) in a brief way. Introducing 
the OV function 

Vopt(x) := V o( Kl + e -( x - Xo )/Sx ~ i + e x /8x 

as well as the effective distance 

A cS x k (t) :=5x\og(±- r e-^^)/ fa t' N "' 
^to Jo 



the s2s-OV model is expressed as 



7 op t(A eff a; fc (t)). 



(Al) 



Let x£ := = n5t) and t>]? := — x%)/5t. Discretization of the effective distance 

A e ff£fc(t) is then expressed as 

A e d ff :4 :=5*log(£ 



n'=0 



n Q + 1 



where n := to/5t. Using the fact that the OV function v opt (x) is written as a limit of a 
function 



ox 

^optO) : = 77 lo g 



(ft 



X -)- g(x-x )/5x /]__)_ e (x— xo— fo<5*)/<5a 



X _)_ g-x /<5x / ]_ _|_ g-(xo+^o<5*)/<5x 



as foptt^) = lim^^o ^optl^X the s2s-OV model (IA1I) in a time-discretized form is given by 



vz = v „ 



: (AixJ), 



(A2) 



which is equivalent to 



;r 



n+l „n 



< + 5x< lot 



-log 



»() 



1 + 

n'=0 
no 

1 + 



-(A**"" -Xo)/&D N -1 



n + 1 



-log(l + e- xo/fa ) 



'"I e -(Ax™-" -x -uo<5t)/5x N _i' 



n'=0 



n + 1 



log(l + 



-(xo+^o^O/^x 



) ■ 



It is straightforward to confirm that the continuum limit 5t — > of the above discrete s2s-OV 
model ( ]A2jl reduces to eq. flAlj) or equivalently the s2s-OV model ([I]). 

Ultradiscretization js| is a scheme for getting a piecewise-linear equation from a difference 
equation via the limit formula 

lim 6x \og(e A/Sx + e B/Sx + •••)= max(A B, ■ ■ ■ ). 

5x^+0 
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In the ultradiscrete limit, the OV function for the discrete s2s-OV model v^ pt (x) goes to 
si™ Q v opt( x ) =: v o P t( x ) = max((D, X ^ J - max^O, % ^ - v J, 

which is nothing but the OV function for the uOV model [8:]. The effective distance A^ s x k 
on the other hand is also ultradiscretized in the same manner: 

A\x n k := lim A* s x n k = -max(-Ax"- n ') = min(AxIT n ')- 

6x— >+0 n'=0 n'=0 

Thus we obtain an ultradiscrete equation 

< = < P t(A^), (A3) 

which is equivalent to 

x k +1 = x k+ maxfo, unn(Ax k ~ n ') - x j - maxfo, min(Ai^ n ') - x — v^bt^ , 

as the ultradiscrete limit of the discrete s2s-OV model (1A2D . We name it the ultradiscrete 
s2s-OV model. 

Now let us see how the s2s-OVCA comes out from the ultradiscrete s2s-OV model (IA3I) . 
Assume that the headway between the cars k and k + 1, Ax k := Ax k — xq, for any k must 
be non-negative, Ax£ > 0, which prohibits car-crash from now on. Then the ultradiscrete 
s2s-OV model (IA3[) reduces to the s2s-OVCA (T5]). Note that the constants xq and voSt 
can be arbitrary in the above derivation. We set Xq = 1 and vo8t at an integer so that the 
positions of the vehicles x k become integers at any time. 
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